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摘要 :本 文 给 出 了 四 种 分 析 类 星体 长 周期 光 变 的 方法 ,用 一 个 模拟 的 周期 信号 y=sin0 检验 这 

四 种 分 析 方 法 ， 结 果 表 明 : 四 、 天 体 光 变 采 样 的 数据 点 个 数 相 对 少 到 一 定 值 时 ，Jurkevich 

方法 ， 时 间 补 偿 离散 传 里 叶 变 换 分 析 方法 (DCDFT) ， 离 散 相 关 分 析 方 法 (DCF) 和 功率 

谱 密 度 分 析 方 法 (PSD) 的 分 析 结 果 不 一 样 ， 获 取 最 短 的 连续 数据 点 后 ，Jurkevich 方法 分 析 

结果 在 四 种 方法 中 可 能 最 为 精确 可 靠 ， 且 计算 方法 简捷 实用 。@、 获 得 了 Jurkevich PMA 

法 的 最 佳 参 数 , 当 m=9 时 分 析 结 果 最 佳 ,6@)、 用 m=9 时 的 Jurkevich 方法 分 析 了 类 星体 3C 279 

KIRA 3C 454.3 天 体 的 光 变 周期 ， 得 出 3C 279 天 体 的 可 能 光 变 周期 为 (2.81+0.54)yr , 3C 

454.3 天 体 的 可 能 光 变 周期 为 457d。 
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长 周期 光 变 是 研究 耀 变 体 性 质 的 重要 方法 之 一 , 它 的 确定 关系 到 耀 变 体 的 结构 和 辐射 等 
相关 的 多 个 物理 量 的 估算 ， 比 如 辐射 区 域 的 半径 、 喷 流 的 多 普 勒 因子 和 黑洞 质量 等 ,长 周期 
光 变 为 理论 模型 的 建立 提供 一 定 的 参数 。 光 变 可 以 被 分 为 长 时 标 光 变 、 中 等 时 标 光 变 、 短 时 
标 光 变 申 。 对 于 长 时 标 光 变 和 中 等 时 标 光 变 来 说 ， 由 于 受到 观测 仪器 、 月 相 和 天 气 等 各 种 原 
因 的 影响 ,很 难得 到 比较 完整 的 光 变 的 观测 数据 序列 所。 在 研究 光学 剧变 类 星体 光 变 周期 性 
时 ， 通 过 研究 长 周期 光 变 分 析 方 法 ， 可 以 获取 周期 光 变 分 析 的 最 佳 参 数 ， 用 有 限 的 实际 观测 
数据 序列 获得 最 佳 周 期 估算 值 。 
通常 分 析 周 期 的 方法 有 : 功率 谱 法 、Jurkevich 方 法 、 小 波 、 结 构 函 数 等 ， 这 些 方法 的 特 
点 是 可 以 对 高 于 奈 奎 斯 特 采 样 定律 均匀 采样 的 数据 进行 可 靠 分 析 B5。 但 是 实际 的 数据 分 析 ， 
特别 是 针对 天 体 观 测 所 获得 耀 变 体 的 长 时 标 光 变 周 期 分 析 中 , 这些 方法 的 使 用 就 受到 许多 条 
件 的 限制 ， 比 如 说 Fourier 分 析 外 要求 连 续 的 等 间隔 采样 ， 缺 失 数据 点 的 处 理 就 引进 了 一 些 不 
真实 的 信息 。 所 以 这 些 方法 应 用 在 天 体 的 光 变 周 期 分 析 中 ， 增 大 了 确定 周期 的 误差 。 而 1 
Jurkevich 等 人 提出 的 Jurkevich 方 法 ， 是 一 种 建立 在 期 待 值 均 方 误差 基础 上 的 统计 方法 ， 对 于 
处 理 不 等 间隔 观测 数据 有 很 大 帮助 中 。 本 文 将 利用 四 种 方法 对 这 一 问题 进行 讨论 ， 并 提出 相 
应 的 改进 , 使 得 在 不 等 间隔 时 间 序 列 的 类 星体 光 变 曲线 中 , 寻找 周期 性 的 分 析 计 算 变 得 更 为 


1 ”类 星体 长 周期 光 变 的 计算 方法 


1.1 Jurkevich 方法 
Jurkevich 基于 对 天 文 测 量 中 的 非 均 匀 测 量 问题 提出 的 一 种 统计 法 四， 曾经 文献 [3] 用 该 
方法 分 析 蝎 虎 天 体 PKS 0735+178 的 光 变 特性 。 假 设 我 们 的 观测 样本 数据 为 N 个 ，Xi 为 单 


YEE, X 为 所 有 测量 值 的 平均 值 ， 严 为 测量 数据 样本 方差 ，S 为 样本 标准 偏差 ， 则 有 ; 
pm 1 N — 
Te da, (00 pg-XxLnNy 0 S-JV?^/(N-1). (3) 
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若 样本 划分 为 m 组 ， 则 第 /组 的 统计 参数 应 为 : 


= ] 7 m, —— 
X,-—XX,. (0 Vi-Xxi-mXj. © S?=V?/(m-1). © 
m, i=l i=l 


对 应 m 组 的 总 方差 为 : 
Vass, . (7) 


显然 ，V < 斑 。 上 述 公式 实际 上 表明 ， 当 自 变 量 相同 的 ?个 样本 计算 后 则 变 为 一 个 样 
本 。 如 果 分 组 的 平均 值 相 等 ,那么 最 后 的 期 待 值 方差 将 等 于 分 组 的 平均 方差 。 如 果 分 组 的 平 
均值 不 同 ， 那 么 最 后 的 期 待 值 方差 将 超过 分 组 的 平均 方差 四。 

在 周期 分 析 中 我 们 按照 不 同 的 试验 周期 P 把 数据 折 准 ， 然 后 把 折 闭 好 的 资料 按 其 相位 从 
小 到 大 排列 ， 然 后 分 成 m 组 。 若 给 定 一 个 观测 时 间 为 :f， 试 验 周期 为 P， 则 其 对 应 的 相位 定义 
JU . 


d(t) = mod 
EHARA), t ARMM. EI BV? 04057725 
V2, ， 如 果 所 得 的 实验 周期 等 于 实际 周期 ， 则 歼 的 值 将 会 达到 最 小 ， 由 试验 周期 已 SV? 的 
关系 曲线 中 的 成 最 小 值 对 应 的 时 间 即 可 得 到 样本 的 周期 。 

1.2 时 间 补 偿 离 散 傅 里 叶 变 换 分 析 方 法 
时 间 补 偿 离散 传 里 叶 变换 方法 是 计算 光 变 周期 最 常用 的 方法 之 一 ， 此 方法 最 早 是 由 
Ferraz-mello 在 1981 年 提出 的 外 ,在 过 去 的 研究 中 文献 [10] 用 该 方法 分 析 PKS 1510-089 红外 
光 变 周期 。 通 过 对 ] 、sin wt 、coOs@t 作 Gram-Schmidt 正 交 化 ， 得 到 3 个 正 交 向 量 ， 将 数 

据 投 影 到 3 个 正 交 向 量 上 就 得 到 了 频谱 。 具 体 过 程 如 下 : 


-t 
1) (8) 
P 


H,=1; H,=cos@t H, =sin æt (9) 
正 交 化 后 : (10) h, Sa. 
h =a H, -aa (hs, H) C11) 
h, =a,H, —a,a,(h),H,)—a,a,(h,,H,) (12) 


小 括号 表示 求 两 个 向 量 的 内 积 。 然 后 根据 以 下 关系 确定 ho. his M: 
(hM) = (h.h) E (h,,h,) -] 


(13) 


在 不 均匀 采样 的 情况 下 ， 加 权 的 时 间 补 偿 离散 傅立叶 变换 (DCDFT): 
F(o) = Cf, h, + ipd, 2 
在 许多 类 星体 的 观测 中 ， 观 测 数据 f (6) 的 精度 各 不 相同 ， 考 虑 到 这 个 问题 ， 引 进 权重 
方程 w = wu) 重新 定义 内 积 : 


(2,,2,) = 2 0,81(t;)8,(t;) (15) 
在 内 积 中 引入 权重 后 得 到 : 


ajo = Xa, 


a; =Z o,(cos x)? - a, (La, cos x, )? 17) 
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a5. =} o, (sin x;) -aê (X@; sin x;)? -a [X 9; sin x; cos x; (18) 
-a (Xo; sin x; (Y. cos x;)]' 
在 内 积 中 引入 权重 后 得 到 回归 系数 : 
€ 70 (19) 
c, — 042,0, f (t;)cos x, (20) 


c, =4,} 0, f (t,)sinx, - a.a;a; [2 0; sin x, cosx; — a; (3:6, sin x, (Z @; cosx;)] (21) 


频率 co 处 的 强度 由 下 式 给 出 : M 
I(@)=cq, +c; (22) 
由 线性 回归 理论 可 知 0< T(@)<O, RF: 
Q-(f.f) -Xof(t) (23) 
利用 这 一 性 质 ， 引 进 标准 化 因子 : 统计 量 S(@®) = 7(O)/O ， 我 们 称 这 一 量 为 谱 相 关系 
数 ， 对 于 所 有 频率 办，0<S(oO)<S1 
1.3 离散 相关 分 析 方 法 
离散 相关 分 析 方 法 (DCF) 由 Edelson 和 Krolik (1988 年 ) 用 来 分 析 两 组 离散 数据 相关 
性 的 方法 之 一 , 文献 [11] 用 此 方法 分 析 PKS 2155-304 的 光 变 。 该 方法 可 以 指出 两 个 变量 时 间 
序列 的 相关 性 和 时 间 延 迟 ， 并 可 应 用 于 周期 性 分 析 0]。 该 方法 的 具体 的 步骤 如 下 : 
首先 计算 两 组 数据 的 离散 相关 函数 值 ， 如 数组 a, 和 bp ， 则 离散 相关 函数 值 为 : 
(a; -a)e, -3) 
oe (24) 
a Lb 分 别 为 两 组 数据 的 平均 值 ，c。 和 oj 分 别 为 对 应 的 标准 偏差 。 
其 次 计算 DCF(C) 值 。 通 过 时 间 延 迟 Am =r — 1, 把 两 组 数 联 系 起 来 ， 假 如 时 间 延 迟 为 rz， 在 区 
间 r+tAr/2 中 有 M 个 At;,， 则 DCF(D) 值 为 


DCF(r xr LUD CFy (2) (25) 
再 次 离散 相关 函数 的 误差 为 
o(c)= = (luper, -pcr(cp Ys (26) 

对 于 所 得 到 的 离散 相关 图 , 如 果 峰 值 在 零 的 右面 , 表明 数组 a, 早 于 数组 b, 的 变化 ,反之 ， 
数组 a, 就 会 迟 于 数组 b。 
14 ”功率 谱 密 度 分 析 方 法 

功率 谱 密 度 的 定义 [9 是 如 果 u(1) 是 一 个 可 以 进行 伟 里 叶 变 换 的 函数 ， 则 

u(v)=Su(te2" (27) 
因为 u(t ) 是 实 函 数 ，u(v ) 一 般 来 说 是 一 个 复 函 数 ， 它 们 之 间 满 足 Parseval 公式 


UDCF, = 


xw'()-X zhC)| (28) 
车 ult) 表 示 光 谱 ， 则 等 式 左 端 表示 u(t) 在 (-o%0，oo) 上 的 总 能 量 ， 上 式 右 端 积分 中 的 被 
积 函数 被 称 之 为 能 谱 密度 ， 它 是 一 个 非 负 实数 ， 表 示 了 单位 频率 所 具有 的 能 量 。 在 实际 情况 
下 从 数学 意义 上 讲 大 多 数 w({) 是 不 能 进行 信里 叶 变 换 的 ， 如 果 使 用 堆 尾 本数-(!) 来 截取 
u(t), 
t) <T 
a-lo » 
0, >T 
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那么 对 于 持续 时 间 有 限 的 截 尾 函数 ur() 而 言 是 可 以 进行 傅 里 叶 变 换 的 ， 变 换 式 : 
uiv)» Zulte” = Ful(t)ere (30) 


同样 它 也 满足 Parseval 公式 : 
Yu3(t)= Xl (vf (31) 
联系 截 尾 函 数 ， 将 上 式 两 端 除 以 27 了 并 令 7T 一 % ， 得 到 
1 
Srl) -2-5 |ur( (v (32) 
与 能 谱 密度 的 定义 相对 应 ， 上 式 右 端 被 积 函 数 称 为 功率 谱 密 度 ， 记 为 : 
2 
i) 
P(v)-ii 
(v)= lim 77 


从 表达 式 中 可 见 功率 谱 密 度 是 一 个 非 负 实数 ， 从 整个 功率 谱 密度 的 推导 可 以 看 出 ,功率 
谱 密 度 是 一 个 频 域 中 的 量 ， 它 直接 反应 了 在 频 域 中 不 同 频率 所 对 应 的 值 。 


(33) 


2 ”天文 观测 数据 中 周期 信号 的 模拟 检验 


2.1 天 文 观测 周期 信号 的 模拟 检验 结果 

为 检验 上 述 四 种 研究 方法 的 可 靠 性 , 我 们 用 一 个 模拟 的 周期 信号 作为 天 文 观 测 数 据 来 分 
析 上 述 四 种 方法 。 在 这 里 我 们 以 27x 为 周期 的 正 弱 函数 = sin 0 , 来 检验 分 析 四 种 研究 方法 
的 准确 性 ， 所 使 用 的 单位 为 弧度 。 在 实验 中 ， 我 们 选择 0rad 为 起 点 ， 步 长 为 0.1rad ， 不 
同 的 数据 点 个 数 共 取 了 15 组 ， 所 有 研究 方法 考虑 噪声 等 因素 的 影响 。 
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| 1 Jurkevich 方 法 对 sin 函 数 不 同 数据 点 的 周期 性 分 析 


Fig.l Periodic analysis of sin functions’ different data points by using Jurkevich method 


图 1 为 检验 Jurkevich 方 法 分 析 天 体 周 期 的 可 靠 性 ， 在 这 里 我 们 取 80-360 个 数据 点 ， 每 组 
增加 20 个 数据 点 ， 通 过 Jurkevich 方 法 分 析 得 到 不 同 15 组 数据 点 的 玉 一 PP 曲线 ， 从 而 得 到 的 
分 析 结 果 为 (6.30+0.02)rad 。 
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2 ”时 间 补 偿 离 散 傅 里 叶 变换 分 析 方 法 (DCDFT ) 对 sin 函 数 不 同 数据 点 的 周期 性 分 析 


Fig.2 Periodic analysis of sin functions' different data points by using DCDFT method 


nm 


图 2 为 检验 时 间 补 偿 离散 傅 里 时 变换 分 析 方 法 (DCDFT ) 分 析 天 体 周 期 的 可 靠 性 ， 在 这 
里 我 们 取 80-360 个 数据 点 ， 每 组 增加 20 个 数据 点 ， 用 时 间 补 偿 离 散 傅 里 叶 变 换 分 析 方 法 
(DCDFT ) 分 析 得 到 不 同 15 组 数据 点 的 DCDFT-Frequency 曲 线 ， 从 而 得 到 的 分 析 结 果 为 
(6.33 +0.13)rad ; 
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图 3 ”离散 相关 分 析 方法 (DCF) 对 sin 函 数 不 同 数据 点 的 周期 性 分 析 


Fig.3 Periodic analysis of sin functions' different data points by using DCF method 
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3 为 检验 离散 相关 分 析 方 法 (DCF) 分 析 天 体 周期 的 可 靠 性 ， 在 这 里 我 们 取 80-360 个 
每 组 增加 20 个 数据 点 ， 用 离散 相关 分 析 方 法 (DCF) 分 析 得 到 不 同 15 组 数据 点 的 
elay 曲 线 ， 从 而 得 到 的 分 析 结 果 为 (6.287 土 0.088)rad ; 
/ \ 
| \ \ | | 
| \ "n | mE mu \ 


4 ”功率 谱 密度 分 析 方 法 PSD)〉 对 sin 函 数 不 同 数据 点 的 周期 性 分 析 
Fig.4 Periodic analysis of sin functions’ different data points by using PSD method 


4 为 检验 功率 谱 密 度 分 析 方 法 (PSD) 分 析 天 体 周 期 的 可 靠 性 ， 在 这 里 我 们 取 80-360 


个 数据 点 ， 每 组 增加 20 个 数据 点 ， 用 功率 谱 密 度 分 析 方 法 (PSD) 分 析 得 到 不 同 15 旨 
的 Power-Frequency 曲 线 ， 从 而 得 到 的 分 析 结 果 为 (6.287 土 0.045jrad 。 


图 5 
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日 数据 点 


sin 函 数 用 时 间 补 偿 离 散 传 里 叶 变换 分 析 (DCDFT) 方法 、 离 散 相关 分 析 (DCF) 方 法 、Jurkevich 分 析 方 法 和 功率 谱 密 


度 分 析 PSD) 方法 进行 周期 分 析 时 要 求 的 最 低 数据 点 


Fig.5 The minimum data points required for the sin function to perform periodic analysis using the DCDFT method, DCF method, 


Jurkevichmethod and PSD method 
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由 图 5 分 析 得 到 Jurkevich 方 法 、DCDFT 方 法 、DCF 方 法 和 PSD 方 法 周期 需要 的 数据 点 最 


(KAY 5273504 120, 10047-16047. 3] 
2.2 ”利用 天 文 模拟 数据 寻找 Jurkevich 方 法 的 最 佳 


取 最 短 的 连续 数据 采样 后 ，Jurkevich 方 法 最 为 有 效 。 


在 利用 Jurkevich 方 法 分 析 周 期 时 ， 分 别 取 360 个 、720 个 和 1500 个 数据 点 时 ， 同 时 取 分 组 


数 1-12 组 ， 取 试验 周期 P= 2x，, 令 步 长 为 0.1lraq . H 
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Fig.6 The y — P curve of sin, The number of groupings on the left are 1 to 6 groups,and the grouping on the right are 7 to 12 groups 


如 图 6 所 示 ， 当 取 360 个 数据 点 时 ， 分 纪 
HBA > OM, K 
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第 11 组 和 第 12 组 可 能 是 出 现 了 倍 周期 , 即 可 能 
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图 7 


sin 函数 的 v2 pc 


分 组 数 为 1~6 组 , 右 图 


日数 从 1 组 逐渐 分 到 6 组 的 情况 下 ， 要 求 不 等 间隔 的 
此 周期 性 可 忽略 ; 当 周 期 从 7 组 至 12 组 的 时 候 ， 在 第 10 组 、 
刚好 是 第 一 个 周期 的 重复 出 现 。 即 使 出 现 极 小 


aD Ba TTT TTT 


BOO SOROS PER PSO SOS PP PP 
CNAMRWOONEMBOONAMOOCOMAMACONAMAOOCONAAWDS 


DTTUTUTUIUTUUTUUqgT 


Fig.7 The V2 — P curve of sin, The number of groupings on the left are 1 to 6 groups,and the grouping on the right are 7 to 12 groups 


当 取 720 的 数据 点 时 ， 如 图 7 所 示 ， 分 组 数 从 第 1 组 逐渐 分 到 第 6 组 的 情况 下 ， 同 样 没有 周 
期 性 ;， 当 周期 从 7 组 至 12 组 的 时 候 ， 在 第 10 组 和 第 12 组 出 现 了 倍 周期 ， 即 可 能 刚好 是 第 一 个 


周期 的 重复 出 现 ; 周期 在 第 11 组 出 现 了 明显 的 伪 周 期 ， 因 此 第 9 组 周期 性 的 分 析 结 果 最 佳 。 
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图 8 sin 函数 的 ve P 左 图 分 组 数 为 1~6 组 , 右 图 分 组 数 为 7~12 组 


Fig.8 The J’? — P curve of sin, The number of groupings on the left are 1 to 6 groups, and the grouping on the right are 7 to 12 groups 


如 图 8 所 示 、 当 取 1500 个 数据 点 时 ， 分 组 数 从 第 1 组 逐渐 分 到 第 6 组 的 情况 下 ， 周 期 性 不 
明显 ， 所 以 这 几 组 的 周期 可 忽略 研究 。 当 周期 从 第 7 组 至 第 12 组 时 ， 在 第 8 组 、 第 10 组 和 第 12 
组 出 现 了 倍 周期 ， 第 11 组 的 周期 性 显然 也 没有 第 9 组 明显 。 总 的 来 说 ， 周 期 性 的 分 析 结 果 最 
佳 的 应 该 为 第 9 组 。 以 下 将 用 两 个 多 观测 数据 的 源 来 验证 此 结果 。 


3 ”计算 并 分 析 3C 279 和 3C 454. 3 的 光 变 周期 


3.1 3C 279 的 数据 点 及 光 变 周期 
本 文 研 究 的 3C 279 B 波 段 和 3C 454.3 B 波 段 数 据点 主要 是 从 (http: //www.astro.yale.edu/ 
smarts/glast/home.php) 网 站 获取 。 


从 2008 年 到 近期 ， 由 以 上 网 站 可 搜寻 到 类 星体 3C 279B 波 段 的 数据 点 就 有 661 个 ，B 波 段 
的 光 变 曲线 如 图 9 所 示 ， 可 以 看 出 有 几 次 大 的 爆发 ， 最 亮 时 mp=14.9mag， 光 变 曲 线 中 有 6 个 
以 上 周期 存在 ， 这 些 均 满 足 Jurkevich 方 法 中 确认 长 周期 存在 的 必要 条 件 03]。 
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do 3C 279 在 B 波 段 的 光 变 周期 
Fig.9 The light curve of quasar 3C 279 in B band 
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410 ”Jurkevich 方 法 分 析 3C 279B 波 段 的 光 变 周期 
Fig.10 Jurkevich methed analysis of the variability period of quasar 3C 279 in B band 


利用 Jurkevich 方 法 分 组 数 m=9， 图 10 可 得 到 3C 279B 波 段 的 光 变 周 期 (2.81+0.54)yr 。 周 
NI P. 5 V2 的 关系 曲线 中 有 很 多 的 极 小 值 ， 为 了 有 效 地 判断 周期 己 5 V2 的 关系 曲线 中 所 得 
出 周期 的 真实 性 ， 参 考 文献 0 给 出 了 较 好 的 判 据 : 

ub V? 是 进行 归 一 化 的 值 ， 在 归 一 化 的 3,-P, 图 中 ,车 WV? =1.0, ， 则 有 伍 90， 即 在 样 
本 数据 中 没有 表现 出 有 周期 性 的 存在 。 若 太 ? = 0, WA, 这 时 样本 中 存在 的 最 大 周期 能 够 
从 图 中 识别 出 来 。 在 通常 情况 下 ， 当 f > 0.5 时 ， 则 样本 数据 表现 出 含有 非常 强 的 周期 性 ， 
当 /<0.25 时 ， 则 表现 为 不 含有 周期 性 。 进 一 步 分 析 y2 图 ， 选 择 曲 线 较 “平滑 ”部 分 
中 的 最 小 深度 和 噪声 做 实验 ， 如 果 “ 平 滑 ” 部 分 的 相对 最 小 爆发 值 比 “平滑 ”部 分 的 
v KEIO. 则 可 以 进一步 讨论 相应 的 断 周 期 性 。 在 计算 中 分 组 数 m 的 数目 十 分 
重要 ， 组 数 m 分 得 越 多 ， 灵 敏 度 越 高 ， 但 每 组 数据 中 有 少数 的 数据 点 会 在 图 中 产 
生 较 大 的 噪声 。 分 组 数 m 选 取 较 多 时 ， 还 将 增加 y2 的 计算 量 。 反之， 有 可 能 寻找 
T 不 了 周期 。 对 3C 27948 V2 - 尸 图 的 分 析 ， 我 们 很 容易 寻找 到 及 及 已 值 ， 并 用 参考 文 


献 09 的 判 据 ， 三 = (1 一 到 JP 进行 检验 ， 结 果 如 表 1 所 示 。 


表 1 3C 279 天 体 周 期 分 析 表 
Tab.1 The table of 3C 279 objects cycle analysis 


2 2 2 
P E f--xyv; TIẸ 
413 0.7006 0.4273 1.13 
829 0.529 0.8604 2.27 
1222 0.3521 1.8401 3.35 


通过 对 上 述 表 1 进行 分 析 ， 我 们 发 现 3C 279 Adm MV = 0.3521, 71.8401, 可 能 的 
光 变 周期 为 (2.81+0.54)yr 。 利 用 Jurkevich 方法 分 组 数 m=9, 我 们 分 析 了 3C 279 的 B 波段 的 
光 变 数据 ， 图 10 显示 了 B 波段 的 分 析 结 果 。 从 表 一 中 可 以 看 出 ,在 满足 判 据 了 的 条 件 下 [44， 
3C 279 在 B 波段 的 P», Ps 与 Pi1 之 间 存 在 着 简单 的 倍数 关系 : P282 Pi, P3=3 Pi 这 说 明 它 们 


之 间 可 能 存在 着 天 文学 倍 频 关系 。 
3.2 3C 454.3 的 数据 点 及 光 变 周 期 

3C 454.3 (PKS 2251 +158, OY091) 是 一 个 比较 亮 的 、 变 化 比较 剧烈 的 类 星体 ， 并 且 在 
光学 和 射电 波段 存在 着 比较 明显 的 相关 性 03。 本 文 的 数据 是 从 2008 年 到 近期 ， 约 765 个 数据 
点 ， 获 得 了 如 图 11 所 示 的 历史 光 变 曲线 。 
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图 11 3C 454.3 在 B 波 段 的 光 变 周期 


Fig.ll The light curve of quasar 3C 454.3 in B band 
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图 12 ”Jurkevich 方 法 分 析 3C 454.3B 波 段 的 光 变 周期 
Fig.12 Jurkevich methed analysis of the variability period of quasar 3C 454.3 in B band 
利用 Jurkevich 方法 分 组 数 m=9, 我 们 分 析 了 3CA54.3 的 B 波段 的 光 变 数据 , 图 12 显示 
了 B 波段 的 分 析 结 果 。 从 图 12 中 可 以 看 出 ，3C454.3 在 B 波段 的 有 三 个 明显 的 极 小 值 ， 这 
意味 在 3C454.3 的 B 波段 光 变 曲线 中 可 能 存在 三 个 可 能 的 周期 ， 它 们 分 别 是 : P,=457d, 
P2=891d 和 PP; 二 1320d。 然 而 ， 根 据 Kidger 等 人 的 结果 ， 要 确定 一 个 周期 ， 数 据 样本 的 时 
间 跨 度 要 超过 周期 的 六 倍 L9。 而 我 们 的 数据 样本 的 时 间 跨 度 大 约 25004， 这 小 于 Po. Pil 
六 倍 ， 故 了，，P; 必 须 排除 ， 并 需要 更 多 的 观测 数据 来 确定 它们 。 但 非常 有 趣 的 是 P，，P 
3 与 Pi 之 间 存 在 着 简单 的 倍数 关系 : P2 Pi P33 P1， 这 说 明 它们 之 间 可 能 存在 着 天 文 
学 倍 频 关 系 。 这 意味 着 3C454.3 在 B 波段 的 光 变 曲线 中 可 能 存在 一 个 书 ,=457d 的 真正 的 光 
变 周 期 。 


同样 ， 对 3C 
ROM: f 


454.3 ARI V2 一 已 图 的 分 析 ， 我 们 很 容易 寻找 到 及 及 已 值 ， 并 用 参考 文 
=(1 -成 )/V? 进行 检验 ， 结 果 如 表 2 所 示 。 


m 


表 2 3C 454.3 天 体 周期 分 析 表 
Tab.2 The table of 3C 454.3 objects cycle analysis 


2 2 2 
P y Fave yin TIE 
457 0.5505 0.817 1.25 
891 0.58 0.724 2.44 
1320 0.4133 1.420 3.67 


通过 对 上 述 表 2 进行 分 析 ， 我 们 发 现 3C 454.3 4 V2 最 小 值 为 0.55305， 广 0.817， 对 应 


的 光 变 周期 可 能 为 P-457d， 这 既 满 足 了 图 12 的 分 析 结 果 ， 也 满足 / 浏 据 0。 


4 ”讨论 与 结论 


本 文 利用 一 个 以 27 为 周期 的 正弦 函数 y = sin 9 来 检验 四 种 研究 方法 的 可 靠 性 ,在 实验 


中 ， 所 使 用 的 单 


位 为 弧度 ， 选 择 0rad 为 起 点 ， 步 长 为 0.1raq ， 不 同 的 数据 点 个 数 共 取 了 


15 组 。 得 到 Jurkevich 方法 的 分 析 结 果 为 (6.3 土 0.017)rad ; 得 到 时 间 补 偿 离散 傅 里 叶 变 换 


分 析 (DCDFT) 


方法 的 分 析 结果 为 (6.331 土 0.130)jrad ; 得 到 离散 相关 分 析 (DCF) 方法 的 


分 析 结 果 为 (6.287+0.088)rad ; 得 到 功率 谱 密 度 (PSD ) 方法 的 分 析 结果 为 


(6.287 土 0.045)rad 。 分 析 结 果 表 明 Jurkevich 方法 、DCDFT 方法 、DCF 方法 和 PSD 方法 


周期 需要 的 数据 


点 最 低 分 别 为 50 个 、120 个 、100 个 和 60 个 。 获 取 最 短 的 连续 数据 采样 后 ， 


Jurkevich 方法 最 为 有 效 。Jurkevich 方法 分 析 结 果 在 本 次 研究 的 四 种 方法 中 最 为 精确 可 靠 ， 


日 此 计算 方法 简 


捷 实 用 


在 实际 应 用 


中 发 现 Jurkevich 方法 十 分 依赖 于 分 组 数 m，m 值 越 大 分 析 结 果 就 越 好 ， 但 


会 产生 比较 大 的 
数据 分 布 偏差 很 


噪声 ， 反 之 则 有 可 能 找 不 到 周期 (1。 如 果 数 据 分 布 不 均匀 ， 则 会 导致 组 间 
大 ， 从 而 影响 获得 实际 的 周期 。 根 据 图 6、 图 7 和 图 8 可 以 看 出 ， 对 于 不 同 


的 m 值 ，V? 与 ,的 关系 曲线 的 倾斜 程度 不 同 ， 周 期 大 的 向 下 倾斜 。 结合 检验 公式 可 知 ， 


m 


户 越 小 了 值 就 越 大 ， 就 更 容易 满足 f 检验 。 由 此 可 知 检验 依赖 于 m 值 的 大 小 。 进 而 利 


用 模拟 数据 寻找 到 Jurkevich 方法 的 分 组 数 为 m=9 组 时 分 析 结 果 为 最 佳 。 最 后 用 分 组 数 为 


m=9 组 时 的 Jurkevich 方法 分 析 了 3C 279 天 体 及 3C 454.3 天 体 的 光 变 周期 ,我 们 得 出 3C 279 
天 体 可 能 的 光 变 周期 为 (2.81+0.54)yr ， 在 文献 [18] 中 分 析 了 3C 279 光 变 可 能 周期 为 
130.641.3d ， 我 们 分 析 得 到 的 可 能 周期 大 约 是 130.6+1.34 的 8 倍 。3C 454.3 天 体 可 能 的 光 变 
周期 为 4574。 在 文献 [19] 中 分 析 了 3C 454.3 光 变 可 能 周期 12.39yr， 大 约 是 我 们 得 到 的 周期 


的 10 Èo 


致谢 : 感谢 编 
这 些 意 见 对 我 们 
了 仔细 修改 。 
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Study on long-period variability analysis methods ofquasars 
YuLian, Zhang Xiongt, ZhangHaojing, LiFuting, XuXiaoling, RenGuowei, WuYuecheng 


(College of Physics and Electronics, Yunnan Normal University, Kunming 650500, China, Email: ynzx@yeah.net) 
Abstract:This paper first introduces four methods for analyzing long-period variability of quasars, 
then useing a simulated periodic signal y = sin 0 to test four methods.The analysis results show: 1. 
When the number of data points sampled by the celestial body is relatively small to a certain 
value,the analysis results of Jurkevich method,DCDFTmethod, DCFmethod andPSDmethodare 
different, After getting the shortest continuous data point,the analysis result of Jurkevich method is 
probably the most accurate and reliable among the four methods, and the calculation method 1s 
simple and practical.2.The best parameters of the Jurkevich analysis method were obtained, The 
best results are when m-9.3.Using the Jurkevich method withm=9 to analysis the variability of 
3C279 is about (2.81: 0.54)years and thevariability of 3C 454.3 is about 457d. 

Key words:Quasars; Jurkevichmethod; DCDFT method; DCF method; PSD 
method;Long-period variability 


